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1 Introduction 

Despite the enormous effort done since longtime, there are only few exactly 
solvable nonlinear field equations in (1 + 1) dimensions. Among these models, 
the Burger equation is a well known example of diffusion nonlinear equation. 

In their remarkable contribution [1], Th. W. Ruijgrok and T. T. Wu pre- 
sented a completely solvable model of a discrete velocity, non-linear Boltzman 
equation (here referred as the RW equation). Basically, the non-linear evolu- 
tion can be linearized via a logaritmic transformation of the Hopf-Cole type. 
In that sense, the RW model has been interpreted to be a generalization of 
the Burger equation [2]. Indeed, the Burger equation describes the evolu- 
tion of the velocity field of a transport process governed by the White noise 
process (WNP), [3], while the RW equation can be viewed similarly as the 
evolution for the coupled velocity fields of a transport process governed by 
the random telegraph process, (RTF) [2]. 

In this paper we present a completely solvable nonlinear reaction-diffusion 
model. The starting point are the Chapman-Kolmogorof equations corre- 
sponding to a dynamical system driven by a sum of the WNP and the RTR 
The nonlinear equations obtained after an inverse Hopf-Cole transformation 
are a couple set of nonlinear reaction-diffusion equations. The exact solution 
of this dynamics is expressible in terms of a linear differential operator acting 
on the logarithm of the convolution product of a Gaussian with the solution 
of the telegrapher's equation. In this way it is possible to construct all the 
solutions corresponding to initial conditions in the physical domain, i.e. such 
that they remain positive, as it might be for a distribution field. The model 
[1] and its solution is restituted when one of the control parameters (i.e. the 
diffusion constant governing the WNP) is set equal to zero. 

2 The model and its solution 

Consider the coupled set of non-linear reaction-diffusion equations : 
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where a, (3, v, /i are positif reals, S = fi + f2 and D = fi — /2. Subscript 
X denotes the derivative with respect to x and the arguments (x, t) of /i and 
/2 have been omitted in the definition of K^ for simphcity. 
Remark . When /i = 0, the model defined by Eqs. (1), (2) and (3) reduces 
to the RW equation [1]. In this case the model has the following properties: 
(1) conservation of the number of particles and energy; (2) nonlinearity; (3) 
positivity of distribution functions and (4) uniqueness of the equilibrium 
distribution for any given density. 

In our case diffusion clearly prevents property (1) and property (4) may 
only be verified with respect to each initial condition. The remaining prop- 
erties are still verified in our case. 

Definition 2.1 We define the following domain of admissible initial condi- 
tions ' 

T{a,(3,fi,v)= {/ix/2GC2(R)xC2(R)|D(x,0) 

(4) 

-i^S,{x,0)-J^S\x,0) + a + (3} 

We denote by G{x, t) the usual density for the WNP, namelly 

and by T{x, t) the corresponding density for the RTP, namelly 

nx, t) = W{x, t)exp{-^t - ^x} (6) 

where 
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where /q and Ji denote modified Bessel functions and 



Wq{x) = exp 



S{x',0)dx' + 
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Wi{x) = \ \d{x, 0) - ^S^{x, 0) - ^S\x, Q) + a + (3\w, 
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Tlie next proposition is tlie main result of this paper. 

Proposition 2.2 The solution of the set of reaction- diffusion equations (1), 
(2) for initial conditions {fi{x, 0), /2(x, 0)) G F, reads 

(10) 
(11) 

(12) 
(13) 



fi{x,t) = 0~\n{N{x,t)) 

f,{x,t) = Onn{N{x,t)) 
with the operators O^ defined by : 
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and 



N{x,t)= / G{x',t)T{x-x',t)dx' 



Proof We follow the idea of the construction presented in [1]. Hence, let 
the operators O^ as defined in Eq.(12). These operators commute. Therefore 
it exists a field F{x,t) such that Eqs. (1) and (2) are equivalent to : 



0+0-F{x,t) = -K^ Lo'F,0+F,-0--^F, 
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with 

-0~F{x,t) = hix,t) (15) 

and 

O^Fix,t) = Ux,t) (16) 

Consider now the field N{x, t) : (x, t) G /? x R'^ -^ R^ defined by : 

N{x, t) = exp{-F{x, t)} (17) 

Introducing Eq. (17) into Eq. (14), after some tedious manipulations, it is 
possible to verify that /i, /2 will satisfy the system (1) - (2) with K^ defined 
as in (3) provided the field N{x, t) obeys the following evolution equation : 

det(-'";°'' _|,3^O-,)a'(-.')=0 (18) 

Eq. (18) is a fourth order hyperbolic partial differential differential equation 
with a probabilistic interpretation, [4, 5]. Consequently, Eq. (18) conserves 
the positivity of its solution. This equation has been first encounterd in [4] 
(see Eq. (21) of this reference) and its solution is given in [5], where it is 
shown that : 

N{x,t) = {G*T){x,t) (19) 

with * denoting the convolution product as in Eq.(13) with : 

|G(x,t)=^2^G(x,t) (20) 
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From Eqs. (10), (11) and (12), we have : 

N{x, 0) = No{x) = exp{-^ J dx (/i(x, 0) + h{x, 0))} (22) 

and 

^iV(a;,0) = iVi(x) = i{/i(x,0)-/2(x,0) + ^^[/i(x,0) + /2(x,0)]}iVo(x). 

(23) 



In order to get G{x,t) as in Eq.(5) it is sufficient to start with the initial 
value in the form : 

G{x,0) = 5{x) (24) 

Now, concerning T{x,t), in order to match Eq. (6.4) of [1] we introduce 
W{x,t) defined by Eq. (6). In view of (19) and (6) N{x,t) will be known as 
far as the solution W{x,t) of the following equation is found : 

^W{x,t)=(^v'^ + aP^W{x,t) (25) 

To integrate (25) one needs initial conditions. The last are related to the 
initial conditions for N{x, t) by 

Wix, 0) = Woix) = exp{ [ ^^ ] x}Noix) (26) 

and hence with the initial conditions for /i and /2. Therefore they may be 
computed using (19)-(24) which leads to (8) and (9). 

Finally, thanks to Eqs.(8), (9) and (26), the solution of Eq. (25) is given in 
[1] and reads as written in Eq. (7). For the initial values (/i(x, 0), f2{x, 0)) G 
r as defined by Eq. (4), Wo{x) and Wi{x) are positive. Hence W{x,t) > 
0,\/{x,t) E Rx i?"*", see Lemma 1.3.2 in [6]. Hence N{x,t) is positive since 
by Eq.(19) it is the convolution of two positive fields. Therefore ln[N{x,t)) 
is well defined. The proof is complete. 



3 Concluding remarks 

In this note, we have shown how to extent the exactly solvable discrete ve- 
locity non-linear Boltzman equation introduced by Th. W. Rujgrok and T. 
T . Wu [1], (the RW equation). In this way we obtain a nonlinear reaction- 
diffusion coupled system of equations for which the exact solutions can be 
computed as far as the corresponding initial conditions are in the positivity 
preserving domain. Our construction is based on the observation that that 
the RW equation and the Burger equation are non-linear evolutions for ve- 
locity fields of linear transport processes, a point of view pioneered in [4] and 
further extended in [2]. Accordingly, the construction presented here can 



be further generalized by considering sets of non-linear equations for veloc- 
ity fields describing a transport process with dynamics governed by a mixed 
noise stochastic differential equations [5]. Therefore, the model presented 
here is only one of the simplest of a whole family to be further investigated. 

From that perspective, the special case sorted out here might appear to 
be an arbitrary choice. However, due to to the simple form of its solutions 
and its reaction-diffusion form, it may deserve special attention. 

In a forthcoming paper we shall study some properties of the solutions of 
this model. 
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